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Short Coherence Length Superconductivity: A Generalization of BCS Theory 
for the Underdoped Cuprates 
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On the basis of the observed short coherence lengths in the cuprates we argue that a BCS-Bose-Einstein condensation (BEC) 
crossover approach is an appropriate starting point for correcting the mean field approach of BCS and, thereby, for addressing 
pseudogap phenomena in these materials. Our version of the BCS-BEC approach is based on a particular Greens' function 
decoupling scheme which should be differentiated from others in the literature, and which yields (i) the Leggett crossover 
ground state (for all coupling constants g, at T = 0) and (ii) BCS theory (for all T <T C over a range of small g). In this paper 
we provide a simple physical picture of the pseudogap phase above and below T c , and review the quantit ative and qualitative 
implications of this theory, which, for the most part have been published in a series of recent papers. cond-mat/0003133 



1. Introduction and Formalism 

While two main alternatives have been suggested 
for addressing the cuprate pseudogap (the phase fluc- 
tuation approach[Jl|] and the "nodal" d-wave quasi- 
particle picture[||]), this paper presents a third alter- 
native which is based on the presumption that one 
should investigate small excursions from BCS with- 
out abandoning it altogether. This philosophy rests 
on the observation that (i) £ is short so that the mean 
field theoretic approach of BCS should not be ex- 
pected to be correct in detail and (ii) on Uemura's ob- 
servation^ that the cuprate superconductors belong 
to a large class of "exotic" materials. While the possi- 
ble role of stripes, spin-charge separation and d-wave 
pairing has received much attention in the cuprate lit- 
erature, here we argue that Uemura's famous plot sug- 
gests that one formulate a theory of superconductivity 
around a more universal approach without invoking 
specific materials-dependent features. 

In this paper we present an extension of BCS 
theory, capable of describing short £ superconduc- 
tors, which is intermediate between a BCS and Bose- 
Einstein condensation (BEC) picture (the "BCS-BEC 
crossover scenario"). At the heart of our approach is 
the distinction between the fermionic excitation gap 
A and superconducting order parameter A sc which 
distinction holds for T ^ 0, both above and below T c . 
The difference parameter (or pseudogap energy scale) 
Ap g — A 2 — A 2 C can be associated with low energy, 
weakly damped, pair excitations of finite momentum. 



It is essential to note that our generalization of 
BCS theory is based on the particular ground state 
of the crossover problem^]: the Leggett state ( = 
rik(wk + u kC k cL k )|0)), which has a BCS-like char- 
acter but should be thought of as applicable to arbi- 
trary values of the attractive coupling constant g. Here 
itk, fk are given by the usual BCS-like expressions in 
terms of the fermionic dispersion E^. The variational 
conditions which ensue from Leggett' s calculations 
are the same as those written below in Eqs. (1) and 
(2), for T = 0. We have shown, based on the work of 
Kadanoff and Martin[|[| (KM), that these self consis- 
tent equations can be derived from a Green's function 
decoupling scheme. We stress here the three impor- 
tant but subtle observations of KM: (i) to obtain BCS- 
like theories, it is sufficient to truncate the system of 
equations so that only one and two particle propaga- 
tors enter and (ii) the correlation functions associated 
with this KM scheme - at the level of the gap equa- 
tions and thermodynamics - are not the same correla- 
tion functions which enter into the electrodynamical 
response, (iii) The pair susceptibility \ which appears 
in the T-matrix, or pair propagator T = g/(l + gx), 
is of the form x = GGo\ this differs from other T- 
matrix schemes (where, frequently, \ — GG), but 
is required to obtain the BCS-like gap equations and 
their thermodynamics [§]. It should also be noted that 
this ground state has 100% condensation for all g, un- 
like that of the non-ideal Bose gas. 

Without any detailed calculations, we can antici- 
pate the form of the generalized BCS theory, which 
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applies below T c . As expected by a natural exten- 
sion of the Leggett ground state variational conditions 
(and, as is consistent with BCS theory) we have 
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Here £k is the bare fermion dispersion, measured 
from the chemical potential \x, — yje^ + A 2 ^^. 
is the quasi-particle dispersion and ip-^ represents the 
symmetry of the pairing state. The new physics of 
the pseudogap phase is embodied in a third equation 
which represents the difference of the two energy gap 
parameters in terms of the number of excited pairs 
(with Bose distribution b(Cl g )), as 



A 2 -AL = A^ =a„VW, 
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where the coefficient of proportionality, do, and the 
pair excitation energy, S! q , can both be determined 
from the above microscopic theory. These calcula- 
tions^^]] indicate that (at small, but non-zero (q, f2)) 
the pair propagator can be approximated by 

T ^ a /[fl - il q + p pair + iT q ] (4) 

for the purposes of calculating the self energy S of 
the fermions (which self energy enters into the re- 
sulting gap equations and thermodynamics). One can 
interpret this theoretical approach as follows. Here 
we go beyond BCS to include self energy £ effects 
associated with the incoherent, finite q pair excita- 
tions. Whereas, BCS is a mean field treatment of the 
particles, the present approach should be thought of 
as a mean field treatment of the pairs. This repre- 
sents the next level in a hierarchy of (superconduct- 
ing) mean field theories, which hierarchy is similar to 
that encountered in magnetic problems. This mean 
field approximation of the pairs is associated with the 
truncation of the equations of motion, so that pair- 
pair interactions are not directly present. Residual 
inter-boson interactions arise only indirectly via the 
self energy of the fermions. As a consequence, both 
above and below T c , the pair dispersion is quadratic 
n q = q 2 /Mpair, as in a quasi-ideal Bose gas. In this 
way, as in the KM paper, the pair correlation function 
is to be distinguished from that associated with the 
collective modes of the order parameter[K[]. 



2. Results 

The results summarized here represent those de- 
duced from solving the coupled equations^] for the 
pair propagator T and its single particle analogue (G), 
along with the fermionic number constraint. 

Pseudogap onset: In our earliest work^ we found 
that the temperature T* at which the Fermi liquid state 
first breaks down corresponds to the onset of a split- 
ting of the single peaked (broadened Fermi liquid- 
like) electronic spectral function into two peaks sepa- 
rated by a (pseudo)gap. This pseudogap state is char- 
acterized by the presence of metastable pairs or "res- 
onances" which effectively reduce the single particle 
density of states. Near, but above T c , the latter takes 
the form of a broadened BCS-like structure. More- 
over, between T c and T*, the wave-vector dependence 
of S (and of the pseudogap) departs from that of a 
strict Lfk symmetry [0| . Note, also, that in the present 
approach we always have T* > T c . 

Superconducting transition: As the temperature 
decreases from T*, the density of meta-stable pairs 
with momentum q continues to increase, until at T = 
T c , a macroscopic fraction with q = undergoes 
Bose condensation. At this temperature, the pair 
propagator or T-matrix diverges as expected from the 
Thouless criterion, so that both p pa ir and T q= o van- 
ish. Moreover, near T c , the inverse lifetime of the q ^ 
pairs, as well as that associated with fermion states 
7 becomes very small[[7|]. This is an important set of 
observations. Once the vicinity of T c is reached, the 
interaction between the long wavelength (soft) bosons 
and the (gapped) fermions becomes weak. This same 
behavior continues to hold below T c so that to leading 
order, lifetime effects can be dropped in Eqs. (l)-(3). 

The dependence of T c on g which results from Eqs 
(l)-(3) is highly non-monotonic, even for the s-wave 
case[0]. At low g, the curve for T c vs g follows 
the BCS result until A pg becomes sufficiently large; 
then T c decreases, with g, thereby reflecting the diffi- 
culty of forming a superconducting state in the pres- 
ence of a fermionic gap. Once n becomes negative, 
T c then increases with increasing g approaching a 
(mass) renormalized, but otherwise ideal BEC limit. 
For the d-wave case, the situation is even more com- 
plex[0. As a result of the extended size of the d- 
wave pairs, the effects of the Pauli principle repulsion 
are more extreme, and T r vanishes well before the 
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system reaches the /i < limit (except in the limit 
of unphysically small densities). In strictly 2d, we 
find T c is always zero. 

Behavior of the superconducting gap equations: 
Equations (l)-(3) were solved to yield the results plot- 
ted in Figure 1, for the case of weak, moderate and 
large g . Also indicated is a "cartoon" characteriz- 
ing the excitations of the condensate at each value of 
g. Above T c , where the computations are more com- 
plicated we used a simple extrapolation procedure^ 
to facilitate the numerics. As can be seen from the 
Figure and from Eq. (3), as the number of excited 
pair states increases, so does the difference between 
the fermionic excitation gap and the order parame- 
ter. These figures represent a natural generalization 
of BCS theory. 

Thermodynamical consequences; mostly C v : The 
pair excitations shown schematically in Figure 1 
necessarily lead to corrections to the quasi-particle- 
derived thermodynamics of BCS theory, via new low 
T power laws. Indeed, it would be difficult to imag- 
ine how fermionic quasi-particles could be relevant to 
the thermodynamics of the BEC limit. That these ex- 
cited pair states have a nearly ideal Bose gas charac- 
ter is a consequence of the BCS-like, Leggett ground 
state; they arise from a mean field treatment of the 
pair propagator. One might imagine an improved ap- 
proach to the strong coupling limit (which includes 
direct pair-pair interactions) would lead to results 
which are more analogous to those of the non-ideal 
Bose gas. It should, however, be noted that the 
cuprates are in the clear fermionic regime where fi 
is essentially Ep, as is consistent with our calcula- 
tions [|l(],[iT]]. Moreover, for the <i-wave case, our anal- 
ysis shows that the BEC limit is virtually inaccessi- 
ble, and concerns about inaccuracies of the KM ap- 
proach, in this limit, appear to be less relevant. What 
is most important is to properly capture the physics 
of the BCS regime, so that small excursions from it 
can be considered in a controlled manner. As a result 
of these pair excitations, in a quasi-2d system such as 
the cuprates we find[[l2|j C v = aT 2 + 7*T. At this 
time, a low T, linear specific heat contribution is well 
documented in the cuprates, although it is not known 
whether it is intrinsic or extrinsic. Elsewhere in this 
journaljotl, we compare our quantitative calculations 
for C v with experiment. 

We have not yet completed a full calculation of the 




Figure 1. Temperature dependence of the gaps. Pair excita- 
tions become progressively more important with increasing 
coupling strength. 

thermodynamical properties such as C v from above to 
below T c . Nevertheless, it can be anticipated on phys- 
ical grounds that C v /T starts to decrease once T* is 
reached and that the overall behavior of that contri- 
bution to C v which is associated with the fermionic 
degrees of freedom, is rather similar to that of a much- 
broadened BCS theory above T c . The density of 
states is depleted precisely at Ep, once the temper- 
ature passes below T c . This long range order-induced 
depletion, is, thereby, reflected in thermodynamical 
properties, such as C v jumps which become progres- 
sively weaker the stronger the coupling g. 

Behavior of the penetration depth: We find[ll 
p"2| ], that the penetration depth at low T, behaves 
as A = A + AT + BT 3 / 2 . Here the second and 
third terms represent respectively the fermionic quasi- 
particle and the bosonic contributions. This additional 
T 3 / 2 bosonic contribution may be difficult to distin- 
guish from the "dirty" d-wave T 2 term, which was 
previously invoked in fits to the data, at the lowest 
T. A quantitative discussion of this term is given 
elsewhere in this journal Jl3|], along with comparisons 
to the data. [The hole concentration (x) dependence 
must be also addressed, as will be summarized be- 
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low]. It should be stressed that this bosonic term is 
responsible for the quasi-universal scaling of the pen- 
etration depth] 11 1, which we have reported. Indeed, 
our calculations indicate that there are systematic de- 
viations from perfect scaling, and the general trends 
for these deviations appear to have been observed by 
the Cambridge group [|T8||. 

The present approach should be contrasted with the 
"d-wave nodal quasiparticle" picture (in its Fermi liq- 
uid rendition [|l4|]) where Landau parameters F\ s are 
thought to be important below T c . In the present con- 
text, these Fig effects are not as primary as is incorpo- 
rating the new (non Fermi-liquid) physics of the pseu- 
dogap state. Ultimately Landau effects can be added 
here, as elsewhere, for detailed fits to data. 

Phase Diagram: In order to incorporate the Mott 
insulator constraint, the band-width or alternatively 
n/m* must be fit to the x-dependence of the zero 
temperature penetration depth flIT|]. In the absence 
of any detailed information about the source of the 
pairing attraction g (which presumably derives from 
Coulomb interactions [|L5[| in some direct or indirect 
form, in the d-wave channel), we assume that g is x- 
independent and fit its ratio to the "bare bandwidth" 
via one adjustable parameter, which is chosen to opti- 
mize agreement with the entire phase diagram [[lO 11 1 . 
With this transcription, one can see that the ratio 
of g to the effective Fermi energy must increase as 
the insulator is approached, so that the bosonic de- 
grees of freedom become more evident with under- 
doping. This prescription provides a quite good fit to 
T*, T c , A(0) over the entire range of x. It also can 
be used as a "tool" for addressing the x and T depen- 
dence of a wide collection of experimental data, in- 
cluding the critical current I c , Knight shift and NMR 
relaxation rate[Q. All of these can be written in 
terms of the contributions from "three fluids": the 
condensate (via A sc ), the fermionic excitations (via 
A), and the pair excitations (via il q , or A pg ). 

Neutron scattering and relation to condensation 
energy: Elsewhere [|l6[|, and in this journal[|l7|| we 
have shown that the incommensurate and commen- 
surate peaks in the neutron cross section can be inter- 
preted as reflecting the d-wave superconducting gap. 
Moreover, above T c , the pseudogap will lead to a 
residue of these peaks, albeit broadened. Within the 
present picture the neutron resonance (commensurate 
structure) and the behavior of the specific heat with its 



discontinuity at and around T c are both consequences 
of the same underlying pseudogap physics, but not di- 
rectly of each other. Moreover, in contrast to BCS 
theory, where one can extrapolate the zero field nor- 
mal state to estimate the condensation energy, here we 
find that because of the non-Fermi liquid nature of the 
normal state (and its instability at T c ) such extrapola- 
tions are problematical. This (along with other fea- 
tures, such as coherent and incoherent contributions 
to ARPES and tunneling, and general finite magnetic 
field effects on the pseudogap) will be discussed in 
more detail in future work. 
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